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ABSTRACT 

Longitudinal vibration of a flexible liquid-filled cylinder was 
studied analytically and experimentally* Experimental observations 
correlated satisfactorily with analysis of the system. Results are 
presented in terms of predicted and measured pressure amplitudes as 
well as the acceleration levels required to produce bubble formations. 
Knowledge of the pressure distribution enables one to estimate the 
vibration limits for which the fluid remains entirely in the liquid state. 
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DEFINITION OF SYMBOLS 


% 


Symbol 

a 



E* 


E" 

f 


G 


Definition 

Fluid acoustic velocity 
Fluid cross-sectional area 
Fluid bulk modulus 
Integration constants 


Longitudinal wave velocity 



Complex wave velocity, 



Cylinder diameter 
Young’s modulus for the container 
Effective longitudinal modulus of elasticity 
Complex modulus of elasticity, E 1 + iE M 
Complex longitudinal modulus of elasticity 
Real component of E* 

Imaginary component of E' 1 ' 


Vibration frequency, 


oo 

2tt 


Net longitudinal force acting on a fluid element 
Acceleration of gravity 

Acceleration amplitude of cylinder, y^oo 2 

8L 

Dimensionless acceleration amplitude, 

g 
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Symbol 



AP 


AP"' 

r 

s 

t 

u 

U 

U m 


DEFINITION OF SYMBOLS (CONT. ) 

Definition 

Threshold acceleration level required to produce bubbles 

Length of the liquid column 

Length of a thin cylindrical fluid element 

Complex amplitude of fluid core motion 

Complex amplitude of boundary layer velocity 

Pressure 

Minimum pressure occurring in the fluid 
Fluid saturation pressure 

Threshold pressure required to produce bubbles 

Ullage pressure 

Oscillating pressure amplitude 

Complex pressure amplitude 

„ . E x 

Ratio, 

B 

Thickness of container wall 
Time 

Longitudinal fluid velocity in boundary layer 

Fluid core velocity 

Maximum core velocity amplitude 
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Symbol 

V 
x 

(i). 

y 

Y 
Y 0 

Yi 

y l 

z 

a 

So 

£ 

E x 

X 

M- 

v 

P 


Definition 

Volume of a cylindrical fluid element 
Depth below surface of fluid 

Location of maximum pressure amplitude 


Radial distance measured normal to cylinder wall 

Longitudinal displacement of fluid cross-section 

Steady component of Y 

Oscillating component of Y 

Amplitude of Yj at cylinder base 

Distance along cylinder wall (see FIG 10) 

Angular deflection of cylinder wall (see FIG 10) 


Oscillating boundary layer thickness 



Circumferential strain of container wall 
Fluid longitudinal strain 

Q 

Wave length, — 


Fluid viscosity 


Fluid kinematic viscosity, 


£ 

P 


Fluid density 
Wall shear stress 
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DEFINTION OF SYMBOLS (CONT. ) 


Symbol 


cr 


4 > 


£2 

a* 

( )*' 


Definition 

Complex amplitude of wall shear stress 
Circumferential stress in cylinder wall 
E u 

Loss factor, rrr 
E 

Angular frequency of vibration 

coL 

Dimensionless angular frequency, -pr~ 
Complex dimensionless angular frequency. 
Complex quantity (time independent) 
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SUMMARY 


Oscillating pressures in a plexiglas cylinder containing 
methanol were measured as it was vibrated parallel to its axis. 
Acceleration levels required to produce bubble formations and the 
locations of these formations were determined. These experimental 
results agree closely with predictions obtained from an analysis of 
the system. However, several discrepancies warrant further 
in ve s tigation . 


INTRODUCTION 


The interaction of fluid dynamics and flexible structures is 
of considerable practical interest in liquid propellant rocket technology. 
Cryogenic propellant storage tanks and lines are known to undergo 
severe longitudinal vibration during flight. Such an environment can 
give rise to two-phase fluid conditions (Ref. 1) that can seriously 
impede the performance of the propellant systems that are required 
to provide rapid fluid transfer (Ref. 2). This study was undertaken to 
develop an under standing of the onset of these two-phase conditions. 


" ^Associate Professor, School of Mechanical Engineering, 

Purdue University, Lafayette, Indiana, NASA-ASEE Slimmer Faculty 
Fellowship Program, June -August 1965. 

**NASA summer student employee, Michigan State University, 
East Lansing, Michigan. 
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Time dependent pressures were predicted from an analysis 
of a flexible liquid-filled cylinder under the influence of longitudinal 
vibration. These predictions were compared with measurements 
obtained using a plexiglas cylinder containing methanol. Knowledge 
of the pressure distribution makes it possible to estimate the limiting 
conditions of vibration for which the fluid remains in the liquid state. 

Beyond these threshold conditions a two-phase fluid state is reached 
in which bubbles are formed. Results of the analysis were used to 
predict the vibration conditions required to produce bubbles and the 
locations of initial bubble formation. These predictions were 
compared with experimental measurements and observations. 

This report is mainly concerned with dynamic conditions 
dominated by single-phase fluid behavior, and the results appear to 
be valid for vibration levels up to those required to bring about the 
onset of bubble formation. 


SURVEY OF LITERATURE 


A sketch of the system investigated is shown in FIG 1. 
Configurations such as this were previously studied by Bleich (Ref. 3). 
Bleich' s investigation is fairly general, because two-dimensional 
effects are considered, and the effect of liquid compres siblity is 
included in the solution of one of the problems discussed. Ponder, 
Blount, and Fritz (Ref. 1) experimentally investigated the bubble 
formations induced by longitudinal vibration of plexiglas cylinders 
(FIG 2). The vibration level required to produce bubbles was deter- 
mined as a function of frequency for several containers and for various 
liquid heights. The experimental and analytical predictions given in 
this report are in good agreement with the observations reported by 
Ponder et al. 

Bleich (Ref. 4), Baird (Ref. 5) and Buchanan, Jameson, and 
Oedjoe (Ref. 6) studied the migration of gas bubbles in vibrating liquid 
columns. This phenomenon is of great interest in the study of two- 
phase conditions. 

The references cited here represent those which are most 
pertinent to the investigation. A comprehensive discussion of all the 
available literature related to this work would, of course, include 
many additional references. 
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SIMPLIFIED ANALYTICAL MODEL 




For ease of physical interpretation and comparison of 
predicted and experimental results, a one-dimensional model was 
analyzed. The following assumptions were imposed: 

1. The fluid behavior is described in terms of longitudinal 
motion, 

2. The fluid remains in the liquid phase with uniform prop- 
erties throughout the volume, 

3. Both the fluid and the flexible cylindrical container remain 
in the elastic range under the loads encountered. 

4. Dissipative effects in the fluid and in the cylindrical wall 
are absent. 1 

5. The liquid remains in contact with the cylindrical wall at 
all times. 

6. The deformation of each cross-section of the cylindrical 
wall is purely radial, and only the hoop stress generated 
by the local radial deformation is significant. 

7. The cylindrical wall is thin, and its inertial effect is 
negligible. 

8. Time dependent pressures in the liquid are not transmitted 
into the gas in the ullage space above the liquid. 2 

9. Inception of bubble formation occurs at a point in the liquid 
when the minimum instantaneous pressure occurring at 
that location decreases to a threshold value, Pq> 


1 The effect of damping is illustrated in Appendix B. 

2 This assumption is easily justified on the basis of the large 
difference in the characteristic impedances (Ref. 7) of the liquid and 
gas phases, primarily because the density of the gas is very small in 
comparison with that of the liquid. 


3 


The equation of motion for the fluid (see Appendix A) is 


* 

v 



d z Y 

ax 2 


(i) 


where Y(x, t) is the longitudinal displacement of a circular cross- 
sectional fluid layer at location x from its equilibrium position under 
zero load conditions* The boundary conditions are 


~E X |“ (0* t) = P U , 


Y(L, t) = Y L e 


icot 


( 2 ) 


The solution of the wave equation (1) is composed of two components; 
a steady component, Yq(x), and an oscillating component, Yi(x, t), so 
that Y(x, t) = Y 0 (x) + Yj(x, t). The equations and boundary conditions 
for Yq and Yi are . 


C 2 


d z Y 0 


g = o 


(3) 


-E x ^ (0) = P u , Y 0 (L) = 0, 


(4) 


and 


9 2 Yi _ ^.2 a 2 Yi 

"aF" “ c 1?" 


(5) 


-E x ^( 0 , t) = o , 1 


Yi(L, 


t) =Y L e 


iojt 


( 6 ) 


Addition of equations (3) and (5) yields the original wave equation (1), 
and addition of the boundary conditions (4) and (6) gives results consis- 
tent with the original boundary conditions (2). 

Adding the solution of equations (3) and (4) to that for equations 
(5) and (6) gives 


* 


T 


This condition is imposed in accordance with assumption 8. 
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P PS 

Y(x, t) = ^ (L - x) + — (L 2 - X 2 ) + 


COS 


Cjl)X ‘ 


Yl 


cos 


ooL 


icot /7X 
e • (7) 


9 Y 

The pressure distribution is obtained from P = -E x — which gives 


P(x, t) = P u + pgx + 



( 8 ) 


where g L is the acceleration amplitude of the base. The steady 
component of pressure at any depth x consists of two parts, the ullage 
pressure, P u , and the hydrostatic pressure, pgx. 

At very low frequencies (— « l), sin — — — and cos — 1. 

\ c / c c c 

For this condition the oscillating component of pressure, like the hydro- 
static pressure, is linearly related to the depth and 

P = P u + p(g + g L e 10jt )x (9) 


In this case the minimum instantaneous pressure at the bottom of the 
container approaches zero when — = fl + ^ j . The corresponding 

minimum instantaneous pressure at other points in the system is given 
by P u (l - )• It can be seen from this example that very low instan- 

V L / gL 

taneous pressures in the system are obtained only if exceeds unity. 

g 

In addition, the lowest instantaneous minimum pressure occurs at the 
bottom of the container. 


At higher frequencies 



the amplitude of the oscillating 


pressure is periodic in x, and the position of the lowest minimum instan- 
taneous pressure moves upward from the cylinder base. For purposes 
of conciseness and generality, equation (8) is rewritten in dimensionless 
form as 
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(10) 


p - p 


U 


PgL 


G 


£2 cos £2 


sin£2 — e*^ 


the lowest minimum instan- 

X TT 


For cases in which the amplitude of the oscillating pressure 
has more than one maximum ^£2 > 

taneous pressure in the system occurs close to the point £2*— = p , 

Q I ^ 

provided that — is somewhat larger than unity. Minimum instan- 

r cos £2 I 

taneous pressures at points farther from the surface and close to 
other maxima of the amplitude term are larger due to the effect of 
the hydrostatic pressure (note the x/L. term in equation (!()))• 


For convenience the amplitude of the time varying pressure, 
AP, is obtained in dimensionless form from equation (10) as 


' AP" 
Pg L 


x 

sin £2 — 
L 

_ G _ 


£2 cos £2 


Resonance occurs for Q = .... and infinite pressure 

amplitudes are predicted by equation (11) for these frequencies. This 
unacceptable result occurs due to assumption 4 and is prevented in the 
real physical system by the presence of friction and possibly by the 
non-linear behavior of the system under large deformations. These 
effects are not accounted for by equation (1), but are discussed in 
Appendix B. 

To initiate bubble formation in a vibrating liquid, the local 
pressure level must be reduced to a sufficiently low value during a 
portion of the vibration cycle. For purposes of simplicity, it is 
assumed that bubbles will form at a particular location if the minimum 
instantaneous pressure falls below some threshold pressure value, 

Pp, (assumption 9). 

For pure liquids it might be expected that Pp would be close 
to the saturation pressure at the liquid temperature. This idea repre- 
sents a considerable oversimplification. In reality the situation is 
far more complicated (Ref. 7) because the threshold pressure required 
to produce bubbles is dependent upon frequency, surface tension, 
dissolved gas properties, and a number of other factors. 
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According to previous discussion in connection with equations 
(8) through (10), the minimum pressure, P m , occurs at the bottom of 
the cylinder and is given by 



( 12 ) 


for £2< -| provided that G > 1. 

- X 

pressure occurs close to Q— 


Similarly, for 12 > ~ the minimum 
= ^ and is approximated by 


Pu ~ Pm 
PgL 


G 

£2 cos £2 


IT 

2 

£2 


(13) 



» 1. 


The condition for inception of bubble formation in the system 
is obtained by setting P m = Px in equations (12) and (13) which now 
are rewritten as 



(14) 


(15) 


Gf represents the threshold acceleration level required to produce 
bubble formation. Equation (14) gives the acceleration level required 
to produce bubble formation at the bottom of the cylinder for G > 1 

when 12 < —, while equation (15) gives the level required to produce 

" X IT TT Q 

bubbles close to — = ^ when £2 > — (provided that -i — r » 1). 

L 2 £2 2 r cos £2 


EXPERIMENTAL DETERMINATION OF OSCILLATING PRESSURES 


A plexiglas cylinder (4-inch diameter by 3 /8-inch thick wall 
by 26 inches long) closed at its lower end was mounted on an electro- 
dynamic shaker. Four pressure transducers were installed in the 
tube wall flush with the inner surface so that the liquid pressures 
could be measured. These transducers were installed at heights of 1, 

6, 9, and 12 inches from the base. A schematic diagram of the experi- 
mental apparatus is given in FIG 3. 

Methanol was introduced into the cylinder, and the system was 
vibrated at levels low enough to prevent any bubble formations. The 
presence or absence of bubbles could be easily viewed due to the trans- 
parency of the plexiglas container. When bubbles formed, the vibration 
was stopped until the bubbles disappeared. The vibration was then 
slowly increased, but was kept at a level low enough so that the bubbles 
did not reform. Thus, all of the measurements and observations 
discussed in this section are associated with vibration of a liquid 
column in the absence of bubbles. 

The pressure signals were amplified and displayed on a four 
channel oscilloscope either simultaneously or one at a time. For 
making pressure amplitude measurements, the latter was more satis- 
factory, because higher gain could be used to expand a single pressure 
trace over a large part of the available range on the oscilloscope screen. 
The results of these measurements are shown in FIG 4a to 4f with the 
analytical predictions from equation (11). Liquid heights of 12, 18, 
and 24 inches were used to give a range of x/L values for the four 
fixed probe positions. 

All of the results show some scatter, a trend that became 
more pronounced as frequency increased. This can be partially attrib- 
uted to the method of plotting the results for comparison with equation 
(11), which does not account for dissipation effects. No doubt other 
effects, not accounted for by the analysis, had their influence in the 
operation of the actual system from which the measurements were 
taken. For example, the one -dimensional longitudinal wave behavior 
(assumption 1) and the absence of wall inertia effects (assumption 7) 
give way to more complicated patterns of motion (Ref. 3 and 8) as 
frequency is increased. 
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4k 


Under conditions of resonance, equation (11) predicts infinite 
pressure amplitudes. For these cases (FIG 4b and 4f), the analytically 
determined mode shapes are given, and the curves are labeled with the 
symbol co to indicate that they are not drawn to scaie. In both cases 
the experimentally determined mode shapes are in agreement with the 
predicted ones. 

For frequencies somewhat removed from resonant conditions 
(FIG 4a, c, d, and e), the measured pressure amplitudes tend to be 
larger than predicted. The deviation is fairly small for the lowest 
frequency illustrated (FIG 4a). 


FIGURE 5 shows the simultaneous display of signals from the 
four pressure transducers as they appeared on the oscilloscope screen 
for the 18 -inch high liquid column. Proceeding from top to bottom, the 
four pressure traces correspond to x/L values of 0. 33, 0. 50, 0. 67, 
and 0. 95. In FIG 5a the frequency is low enough so that all four signals 

are in phase (FIG 4b). FIG 5b corresponds to a higher frequency 

x 

For this situation a pressure node occurs close to — = 0. 67 (FIG 4f). The 

' l-» 

two upper traces are in phase with each other while the third one indicates 
an extremely small pressure amplitude in accordance with its location 
close to the pressure node. The lower trace is 180° out of phase with the 
upper two as would be expected from FIG 4f. 


Since U = u) — = 



it was necessary to obtain the wave 


velocity C for the particular system used in order to evaluate dimen- 
sionless frequency values for various combinations of vibration frequency 
and liquid column height. An estimate of the wave velocity can be obtained 
from the following approximate property values. 


Density of Methanol (Ref. 9)~50 lb/ft 3 

Bulk Modulus of Methanol (Ref. 10) « 1. 86 x 10 5 psi 

Modulus of Elasticity of plexiglas (Ref. 11) % 2 x 10 5 psi 


[b 

Based on these figures the sound velocity for methanol is^— = 4160 fps. 
However, the wave velocity, C, should be lower than this by a factor of 
1 + — — — in accordance with equation (A 10), 

V D E 

value of C = 1250 fps for the experimental system used. 


This leads to a calculated 
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Since the liquid bulk modulus and the modulus of elasticity of 
plexiglas are subject to considerable variation, the calculation given 
above should not be expected to yield a very reliable result in terms of 
the actual system performance. Because of this, C was determined 
experimentally by varying the frequency until the first resonant mode 

was obtained ^£3 = 

of the pressure signals (FIG 4b). The distinct increase in audible 
sound level within a very fine resonant frequency band made it possible 
to tune the system to this condition very easily. With L = 18 M , the 
frequency obtained was very close to 200 cps. Therefore, since 

£2 = 2tt = — , C = 4fL = 1200 fps for the actual system. Thus, the 

previously calculated value of 1250 fps provided a reasonably good 
estimate for C. 

In summary, the simplified analysis predicted the qualitative 
behavior exhibited by the experimental system. In addition, the order 
of magnitude of the calculated pressures is in agreement with the 
experimentally measured values (except for resonant conditions). 
However, a quantitative comparison of predicted and measured pressure 
amplitudes reveals a fairly large percentage of deviations between the 
two for some of the cases studied. 


The mode shape was readily verified by means 


OBSERVED BUBBLE FORMATIONS 


Initiation of bubble formation in a plexiglas cylinder containing 
methanol was investigated. For each case observed, the vibration level 
was slowly increased until bubbles were sighted. The two major item: 
of interest were the location at which initial bubble formations occurred 
and the vibration level required to produce these formations. The 
results are discussed in terms of the simplified analytical model given 
p reviously . 

FIGURES 6a and 6b illustrate system behavior that was typical 
at the lower frequencies. According to the simplified analytical model, 
bubbles would be expected to form initially at the cylinder bottom. This 
did occur, but it was often accompanied by considerable bubble forma- 
tion at the liquid surface when sloshing was present. These bubbles 
migrated rapidly to the bottom after formation. The upper bubble 
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formation region usually would propagate downward, slowly at first 
and then very rapidly, FIGURE 6a shows initial bubble formations occur- 
ring simultaneously at the bottom (since the pressure amplitude was 
largest there - see FIG 4a) and at the top due to sloshing. FIGURE 6b 
shows the rapid growth of these formations that accompanied the 
violent downward movement of the upper bubble formation region. 
Although bubble formation at the bottom is predicted by the simplified 
analytical model, sloshing and bubble formation at the top are not. 

This complicated surface phenomenon, which occurred at low frequen- 
cies in the present investigation, is in agreement with some of the 
observations reported by Bleich and by Buchannan, et al. 


FIGURES 6c, 6d, and 6e show typical initial bubble formations 

at -p- values close to 1. 0, 0. 5 and 0. 33 for Q = tt and —— , respectively. 

These are the locations of maximum pressure amplitude (FIG 4b, 4d and 
4f) and minimum instantaneous pressure, P m , for these frequencies 
according to the simplified analytical model* 


When Q = there are two locations having the same predicted 

x x 

maximum pressure amplitude (FIG 4f), — = 0. 33 and — = 1.0. However, 

L L 

the predicted minimum instantaneous pressure, P m , in the system 
occurs at the upper location, because the hydrostatic pressure is lower 
there than at the bottom. This can be verified from equation (10) which 
contains a hydrostatic component of pressure proportional to depth 
below the surface. Since the predicted minimum system pressure 
x 

occurs close to — = 0. 33, it would be expected that bubbles would form 
L 

at that location first, and not at the bottom. This was confirmed in 
several experiments duplicating these conditions, and the formation 

x 

of one or two small bubbles at — = 0. 33 usually occurred as illustrated 

L 

in FIG 6e. 


On the other hand, the magnitude of the hydrostatic pressure 
was fairly small for the system observed in the present investigation, 

X X 

and the apparent preferential location of — = 0. 33 as opposed to — = 1.0 

did not always dominate. Occasionally initiation of bubble formation 
occurred simultaneously at both locations. Only once did bubbles 
appear at the bottom before any bubble activity was sighted at 
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— = 0. 33; however, a single small bubble formed at — = 0. 33 and 
l-i L 

started to rise to the surface just as the system was photographed 

(FIG 6f). The two bubble formations at the bottom, which actually 

formed first, are clearly visible. The smaller bubble rising toward 

the surface is not as noticeable. 


It was concluded from the many observations made for Q = — - 

M 


that the hydrostatic pressure effect caused the upper pressure ampli- 
tude maximum to be the preferred location for initial formation of 
bubbles. However, this hydrostatic effect was small and did not 
always prevent initial bubble formation at the location of the lower 
pressure amplitude maximum (at the cylinder bottom). 


The photographs in FIG 6a through 6f indicate how the location 
of initial bubble formations is affected by variations in frequency. The 
threshold acceleration level, Gq 1 , was obtained with the same system 
by slowly increasing the vibration level until bubbles were sighted. At 
this point the acceleration level was read from the vibration meter 
(FIG 3). The experimental results obtained are comparable to previously 
reported data (Ref. 1) and are plotted in FIG 7 and 8, For the system 
used in making these observations, the ratio of ullage pressure (atmos- 

p 

pheric) to hydrostatic pressure at the cylinder base was — — = 25. 


FIGURE 7 shows a comparison of the data with the analytical 
prediction for the low frequency range. The analytical result was 
plotted using equations (14) and (15). It was assumed that the threshold 
pressure, P*j, was equal to the fluid saturation pressure, P s , (Ref. 9) 
for the ambient temperature conditions existing at the time the experi- 

Pt Ps P u / Pt\ 

ments were conducted. This gave zr — = ■=— = 0. 2, and — — [ 1 - — — ] = 

Pu Pu PgE V Pu I 

25(1 - 0.2) = 20. The latter value was inserted into equations (13) and 
(14) to generate the plot indicated by the solid line in FIG 7. Both the 
predicted and experimental Gx curves drop sharply in the vicinity of 

the first resonant frequency M) while exhibiting larger values on 


either side of this frequency. The experimental curve does not reach 
a zero threshold value as does the analytical curve because of damping 
in the system. This effect was not accounted for in the simplified 
analytical model (assumption 4) which predicted a zero value for Gq 


12 


at resonance. In view of the method used to generate the predicted 
G-p curve, the agreement of this result with the experimental data is 
considered very favorable in the low frequency range. 

For frequencies somewhat in excess of £2 = the analytical 

prediction gave results considerably larger than the experimentally 
measured Gp values if Pp were assumed equal to P s . This deviation 
had been expected, because it was not realistic to presume that 
threshold pressure is dominated by the fluid saturation pressure for 
all values of frequency, particularly since the liquid used for the 
experiments was not degassed. The bubbles that formed contained 
large amounts of previously dissolved gas. This was apparent 
because the bubbles would rise slowly to the surface due to buoyancy 
following immediate termination of vibration after the bubbles had 
formed. Vapor bubbles would have collapsed instantly. 

Experimental data obtained over a wide frequency range are 
given in FIG 8 . A plot of results obtained from the analytical model 
is also shown. This plot was made using equations (14) and (15) with 

a value of 15 for ( — p- || 1 - 77 ^-]. With = 25, this 

\Pg L )\ P u/ PgL 

Pt 

corresponds to — — =0.4 and a threshold pressure of about 5. 7 psi. 

*u 

The choice of this single threshold pressure level (approximately 
twice the fluid saturation pressure) is somewhat arbitrary since 
P*p in reality is probably frequency dependent as mentioned above. 
However, the use of a fixed Pq> value to obtain a single analytically 
derived curve makes it easier to compare results of the analysis 
with data obtained experimentally . 1 


P T 

1 A smaller value of — — fits the experimental data better up 

it Pu p T 

to fi = - as shown in FIG 7, and a slightly larger value of 75 — would 
2 

TT 

produce more of a favorable comparison for £2 greater than — 
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The similarities in the predicted and experimentally obtained 
Gj curves are apparent in FIG 8. Both undergo distinct decreases at 

TT 3tT 

the two resonant frequencies, Q = — a **d > an d exhibit larger values 

on either side of these two frequencies. Due to damping, the experi- 
mental curve never reaches a zero threshold value. The analytical 
curve, on the other hand, touches the horizontal axis at both resonant 
frequencies. The experimental curve is considerably below the analyt- 
ical curve from £1 = 2. 5 to 4. This is in agreement with the fact that 
measured pressure amplitudes in this range were generally larger 
than those predicted by the analytical model (FIG 4c and 4d). Another 
factor that may explain the deviation in this range is that the actual 
threshold pressure of the fluid may vary considerably with frequency, 
as was previously mentioned. 

The distinct dip in the experimental G-p curve in the vicinity 
of = tt does not agree with the trend predicted from the analysis. 
Possible explanations are: 

(a) The simple one -dimens ional analytical model may not be 
very accurate in this frequency range. 

(b) The damping properties of the plexiglas container may 
have undergone severe variations within this small 
frequency range in a manner which is typical of some 
polymers (Ref. 12). 

(c) The actual threshold pressure of the fluid may vary 
significantly in this frequency range. 

Further effort is needed to determine which of the above possibilities, 
or other contributing factor not listed here, is responsible for the 
behavior of the observed experimental system in the region close to 
' S2 = ir. 


In summary, it was concluded that, except for the surface 
phenomena occurring in the low frequency range, the locations of 
initial bubble formations were predicted reasonably well. The 
predicted threshold acceleration levels produced results corresponding 
to the measured values with respect to order of magnitude as well as 
a significant portion of the general qualitative trends. Very good 
correlation with experimental data was obtained at low frequencies 
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when the analytical prediction was based on a threshold pressure 
equal to the fluid saturation pressure. At higher frequencies, however, 
the quantitative predictions from the analytical model were dependent 
upon the uncertain choice of an appropriate value of threshold pressure, 
and the prediction of pressure amplitudes was not as accurate. Only 
one very distinct trend exhibited by the measured data was not predicted 
by the analytical model. 


CONCLUSIONS AND RECOMMENDATIONS 


In this study many of the experimental observations correlated 
reasonably well with results obtained from an analysis of the system. 
Certain discrepancies remain, however, and to explain them additional 
experimental data as well as more extensive analyses are required. 
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APPENDIX A 


DERIVATION OF EQUATION OF MOTION 


Consider an element of the system as shown in FIG 9. A 
force balance on the cylindrical wall gives 

PD dx = 2crsdx . (Al) 


For an elastic wall, cr = Ee, where e is the circumferential strain of 
the wall. If AD represents the increase in diameter due to pressure, 

P, then e = . Equation (Al) now can be rewritten as 


D 2 P 
" 2sE * 


(A 2) 


Now consider a thin cylindrical slice of fluid of length, i , and volume, 
V =-D i. For small deformations the change in volume of the fluid 
element due to pressure, P, is given by 


AV => 2 / 
4 





(A3) 


Combining equations (A2) and (A3) gives the longitudinal strain of the 
fluid as 


e 


x 




(A4) 


where E x , the effective longitudinal modulus of elasticity of the fluid 
column, is defined as 



(A5) 


Thus, the longitudinal stiffness of the column is dependent upon the 
stiffness, B, of the fluid and the stiffness of the cylindrical container, 
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— E. The ratio of the longitudinal stiffness to the fluid stiffness is 

E x 

less than unity in general and is given by -g- = r =| 


1 + 


B 


(H 


The displacement of a cross-sectional fluid element from 
its equilibrium position under zero load is Y(x, t). The local longi- 

3 Y 

tudinal strain is e x - 7 — . and in accordance with equation (A4), the 

3 Y 

pressure of the fluid is P = - E x — — . The instantaneous velocity of 
8 Y 9X 3 2 Y 

the element is , and its acceleration is — Newton's law is 
written as 


9 2 Y 


pA dx at 2 “ Fx 


(A 6 ) 


where F x is the net force on the element in the x direction. FIGURE 
10 shows the forces acting on a fluid element. Since P may vary with 
x, the diameter of the cylindrical container may vary also, and a 
longitudinal component of force is generated by the tapered cylindrical 
wall. This longitudinal component is PirD dz sin or. Now, sin a = 

1 tb 0 \ /dx\ _ _ 1 _ (bA^ ^dxj 


2\3x/\dzy trDySxy^dz 




becomes P 


V* 


I dx. 


so the tube wall longitudinal component 
The net force in the x direction is 


F x = 


pgA dx + 
(gravity) 


PA 

(force on 
upper surface) 

P 3A 

bx 


PA + 


— J 


(force on 
lower surface) 


= pA dx 


1 3P“ 
g 'p 8 x 


(A7) 


(longitudinal component of 
force from the tube wall) 

3Y 

E x — into equation (A7) and inserting the resulting 


Entering P = 
expression for F x into equation (A 6 ) gives 


9 2 Y 
3 1 2 

(inertia effect) 


aZY 

C + g 
(elastic force) (gravity force) 


(A9) 
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This is the wave equation where 
velocity of the system. 



is the longitudinal wave 


Note that the wave velocity, C, is considerably less than the 


fluid acoustic velocity 




jwhen the fluid stiffness, B, is signifi- 

g 

cantly greater than the wall stiffness, — E. This can be seen by- 
inserting the expression for E x into the relation for C which gives 



APPENDIX B 


FURTHER DETAILS OF SYSTEM BEHAVIOR 


Dissipation of energy can occur in the fluid due to viscosity 
and in the cylindrical wall due to viscoelastic behavior. For conditions 
encountered in this investigation, the fluid viscous boundary layer is 

3 Y 

very thin. 1 The core velocity, U, is equivalent to— ^in terms of 


previously used notation. It is valid to neglect convective terms (of 


the order of U 



in comparison with local accelerations (of the order 


« caU 


m 


where U 



is the maximum core velocity amplitude in the system, 
the above inequality can be rewritten as U m « • 


Since 


The result of the order of magnitude analysis carried out above 
implies that convective momentum effects in the fluid are unimportant if 
the longitudinal wave velocity is considerably greater in magnitude than 
the longitudinal fluid velocity oscillations. This condition was easily 
fulfilled in all of the experiments carried out. In view of this approxi- 
mation the boundary layer (Ref. 13) momentum equation is 


9u _ 9U [ 9 z u 
9t ~ 9t + V dy 2 


(Bl) 


where u is the longitudinal velocity in the boundary layer. 


1 The oscillation boundary layer thickness, or depth of 
penetration (Ref. 13) as it is sometimes referred to, is calculated from 

&o =*]— . For methanol the viscosity at room temperature is about 

\ Cd 

0.55 centipoises. Since the density is 50 lb/ft 3 * (Ref. 9), the kinematic 
viscosity v is 0.74 x 10" 5 ft^/sec. Entering this value into the expres- 

sion for 6 0 gives 6 0 = 0. 001 inch for a frequency of 170 cps. 


19 


Now, the oscillating portion of tne longitudinal core motion is 

t 9 Y j * i t 

written as Y x (x, t) = M v (x)e . Also, U = = icoM^e . In a similar 

manner the boundary layer velocity for a given value of x is written as 
u = N*(y)e 1C °^. Equation (Bl), written in complex notation, is 




The solution is 


N' ,: = iuM* + Ci exp 


(-ryJr) + c * exp ( y 'Jf} 


As y becomes large, must approach iojM'% so C 2 = 0, If the longi- 
tudinal motion of the cylindrical wall follows that of the base, then 

the wall velocity is given by icoYj^ e The fluid velocity at the wall 

is the same, and therefore, io>M* + C\ = uoYl,. With C\ and C 2 
determined, equation (B3) is rewritten as 

N* = icoM * + ia)(YL,- M*) exp ^-y/i-2- ^ (B4) 

where 6 0 is the oscillation boundary layer thickness \/^p • 


The wall shear stress, t, is obtained by differentiating 
equation (B4) and evaluating -p ~ at y = 0. The result, written in 
complex notation, is 


- . ,««(£) 


(Y l~ M*) . 


Equation (A7) is now modified to include the shear force acting on the 
fluid at the wall. The net force becomes 

F x = pA dx^g - + nrDdx. (B6) 

The last term in this equation is obtained from equation (B5). 


To account for viscoelastic behavior, the complex modulus 
of elasticity (Ref. 12) of the cylinder wall is written as E* = E 1 + iE 1 
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! 


For small amounts of damping the real portion is essentially the same 
as the ordinary modulus of elasticity. In this case E* = E(l+ ic|>) where 

E 1 1 

<J> is the loss factor -g-p • The oscillating components of stress and 

strain are obtained from the development given in Appendix A, but with 
E replaced by E* to account for the damping effect. 


Equation (A5) is rewritten as 



(B7) 


The oscillating component of force is obtained from equation (B6) and 
is equated to the inertia term as in equation (A6). With Y x (x, t) 

expressed as M*(x)e 1C °*, the result is 


pA dx(io>) 2 M* = -A dx 


dx 


(WfO 


+ t*ttD dx . 


(B8) 


Inserting t* from equation (B5) gives 


C* 2 d d ^ ~ + to 2 M* + 4(i) 



ca 2 (Y t - M v ) = 0 


(B9) 


=V¥- 


the complex wave velocity of the system. Note 


where C' 1 ' 

1 P 6 

that the last term is negligibly small when « 1, provided that M* 1 ' 

is the same order of magnitude as Yj^. For experiments carried out 

in this investigation, was of the order of 10 3 . It is concluded, 


therefore, that the damping effect due to viscous fluid shear at the 
wall was very small. 1 For this reason the last term in equation (B9) 
is neglected. The solution for the complex pressure amplitude follows 
essentially the same procedure as was used in obtaining equation (11) . 
The result is identical in form and is written as 


1 There is another source of energy dissipation associated with 
viscous action in the core of the fluid due to its periodic compression 
and expansion. In the present study this effect was found to have a 

magnitude even smaller than that of the wall friction. 


i 


I 
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AP* 

Pg L 


X 


(BIO) 


sin J 2 * — 

J-J 

J2* cos Q* 


icot 


In this formulation the oscillating pressure is given by AP" e 
is the dimensionless complex frequency defined as 7 ^. Using the 


n* 


expression given 


in equation (B7), it can be shown that \ f ■ Vy - . 

1 1 + iq> 


The magnitude of the pressure oscillations is obtained from 
equation (BIO) and is written as 


AP* 

Pg L 

G 



n 


1 + irifr 
1 1 + i<(> 


cos 




r<^> 

i<t> 



(B 1 1) 


For the experimental system used in making the pressure measure- 
ments, r was approximately 0. 10. Although the damping properties 
of plexiglas are subject to considerable variation, a loss factor of 
0. 05 appears to be a reasonable value (Ref. 15) at room temperature 
for the frequency range investigated. 

A plot of equation (Bll) is shown in FIG 11 for — = 0. 5 with 

r = 0. 10 and = 0. 05. Predicted pressure amplitudes at the resonant 
frequencies are large, but finite. The limited amount of experimental 
data shown correlates reasonably well with the prediction. In the 
absence of energy dissipation (<|> = 0) equation (Bll) is identical to 
equation ( 11 ), and infinite resonant pressure amplitudes are predicted. 
This unacceptable result illustrates the necessity for including dissi- 
pation effects in any complete analysis of the system. Otherwise, 
realistic predictions in the resonant ranges of frequency cannot be 
obtained. 


To obtain the threshold acceleration levels, it is possible to 
develop relations equivalent to equations (14) and (15) that include 
dissipation effects. It was found, however, that a more straight- 
forward procedure could be carried out using a computer printout of 
equation (Bll). The expression for pressure is first written in a form 
similar to that of equation ( 10 ): 
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g-JP-a -i + 

pgL L 


\pg L ] 


(B12) 


The minimum pressure in the system is equal to P-p at the threshold 
condition. This results in 


gl-ZJjl = 2L _ 

pgL L 


AP*| 

Pg L 


(B13) 


X 


where the — value corresponds to the location of minimum pressure. 
Rearrangement of equation (B13) gives 


^ ■ m . ■»-» *JU a 


AP* | 

Pg L 


(B14) 


The threshold condition corresponds to the position producing the 

lowest value of G from equation (B14). For / ||1 - — 2- ) » 1, this 

)\ P u/ 


occurs very close to the location where the denominator of equation 

(B14) is a maximum. If this position is denoted by(-y-) , then 

W m 




AP* | 
PgB 


(B15) 


maximum 


Inspection of a computer printout of equation (Bll) gave the 
1 AP*| 


maximum values of 


PgL 


and the [-y-) locations. Inserting these 
G V^/m 

values into equation (B15) gave the corresponding threshold accelera- 
tion levels. The result of this prediction is shown in FIG 12 for 
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the same value that was used in plotting the 


result shown in FIG 8. However, in this case energy dissipation is 
accounted for (r = 0. 10, <J> = 0. 05). The predicted values at 
resonance are non-zero in agreement with the behavior of the experi- 
mental data. Despite the improvement, predicted G^ values are still 
low at resonance. Also, the predicted values remain somewhat greater 
than the measured values in the range of £2 = 2. 5 to 4. Thus, based on 
realistic values of the properties of plexiglas, accounting for the 
presence of energy dissipation does not remove the discrepancy in 
this frequency range (FIG 8 and 12). 


In this appendix certain improvements were obtained in the 
predicted results by accounting for viscous energy dissipation. The 
most significant improvements, with respect to the simplified analyt- 
ical model, were obtained at the near resonant frequencies. Certain 
discrepancies between predicted and measured results, somewhat 
removed from the fundamental longitudinal resonances of the system, 
remain to be explained. 
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FIGURE 2. BUBBLE FORMATIONS IN A VIBRATING LIQUID. 

(Reference 1) 
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FIG. 3. EXPERIMENTAL APPARATUS FOR MEASUREMENT OF OSCILLATING PRESSURES 









nalytical: Equation (11) 



FIG. 4c. OSCILLATING PRESSURES WITH O 



Analytical: Equation (11) 



FIG. 4d. OSCILLATING PRESSURES WITH Q 








FIGURE 5 



(a) n=n/2 



lb) A*3n/2 

OSCILLOSCOPE DISPLAY OF OSCILLATING PRESSURES. 
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FIGURE 6a. INITIAL BUBBLE FORMATIONS IN A VIBRATED 
METHANOL COLUMN (£2 = II /4). 
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FIGURE 6b. GROWTH OF BUBBLE FORMATIONS IN A VIBRATED 
METHANOL COLUMN (O = IT /4. 
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FIGURE 6c. 


INITIAL BUBBLE FORMATIONS IN A VIBRATED 
METHANOL COLUMN (S2 = U. 12. 
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FIGURE 6e. INITIAL BUBBLE FORMATIONS IN A VIBRATED 
METHANOL COLUMN ( Q = 3 n/2). 
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FIGURE 6f. INITIAL BUBBLE FORMATIONS IN A VIBRATED 
METHANOL COLUMN (£ 2 = 3 11 / 2 ) 


-Analytical: Equations (14) and (15) 
wlth(P u / pgL)(l-P T /P u ) = 20 

-Experimental: P / pgL = 25 



WITH ANALYTICAL PREDICTION BASED ON P, 




FORMATION 
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FIG. 9 





FORCES ACTING ON A FLUID ELEMENT 





FIG 11 EFFECT OF FREQUENCY ON PRESSURE AMPLITUDE AT x/L = 0.5 
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